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ABSTRACT 

Context. The proximity of core-collapse Supernova 1987 A (SN 1987 A) in the Large Magellanic Cloud (LMC) and its rapid evolution 
make it a unique case study of the development of a young supernova remnant. 
Aims. We aim at resolving the remnant of SN 1987A for the first time in the 3-mm band (at 94 GHz). 

Methods. We observed the source at 3-mm wavelength with a 750-m configuration of the Australia Telescope Compact Array (ATCA). 
We compare the image with a recent 3 -cm image and with archival X-ray images. 

Results. We present a diff'raction-limited image with a resolution of 0^7, revealing the ring structure seen at lower frequencies and 
at other wavebands. The emission peaks in the eastern part of the ring. The 3-mm image bears resemblance to early X-ray images 
(from 1999-2000). We place an upper limit of 1 mJy (2 cr) on any discrete source of emission in the centre (inside of the ring). The 
integrated flux density at 3 mm has doubled over the six years since the previous observations at 3 mm. 

Conclusions. At 3 mm - i.e. within the operational domain of the Atacama Large Millimeter/submillimeter Array (ALMA) - 
SN 1987 A appears to be dominated by synchrotron radiation from the inner rim of the equatorial ring, characterised by moderately- 
weak shocks. There is no clear sign of emission of a diff'erent nature, but the current limits do not rule out such component altogether. 

Key words, circumstellar matter - supemovae: SN 1987A - ISM: supernova remnants - Magellanic Clouds - Radio continuum: ISM 
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1. Introduction 

Supernova 1987 A, located in the Large Magellanic Cloud 
(LMC), was first observed on 1987 February 24. Since then the 
evolution of the supernova remnant (SNR) has been followed 
across the electro-magnetic spectrum, except for the far-infrared 
(FIR) and high radio frequencies. In the radio, since ~ day 1200 
the remnant has been continuously monitored in the 1 .4-9 GHz 
frequency range with the Australia Telescope Compact Array 
(ATCA). At these frequencies, the flux density of the remnant 
has been steadily increasing over time (Manchester et al. 2002), 
in an exponential fashion since ~ day 5000 (Zanardo et al. 2010). 
This reflects the propagation of shocks through the magnetised 
plasma resulting in synchrotron radiation. However, it is as yet 
unclear what is happening at the high-frequency end of the radio 
spectrum, where other emission may become apparent including 
possibly that associated with dust or an emerging pulsar. 

Warm dust (^ 180 K) in an equatorial ring ^ 1''6 in diam- 
eter is linked to mass loss from the progenitor star (Bouchet et 
al. 2006; Dwek et al. 2010). Recently, Matsuura et al. (2011) 
detected a second dust component, of only ^ 18 K, which 
is also seen at 0.87 mm wavelength (Lakicevic et al. 2011). 
Matsuura et al. argued that the cold dust had formed in the ejecta 
of SN 1987 A. Their observations did not have suflftcient angu- 
lar resolution to prove this beyond doubt. While Lakicevic et 
al. (2012) constrained the sub-mm emission to originate from 



within a few arcseconds from the explosion site, they also oflfer 
scenarios in which the dust might not reside in the ejecta. 

The remnant of SN 1987A (hereafter SNR 1987A) was first 
successfully observed at 3 mm wavelength in 2005 with the 
ATCA on short baselines. The remnant was unresolved in the 5'' 
beam; the integrated flux density was 11.2 + 2.0 mJy (Lakicevic 
et al. 201 1). We have now attempted to spatially resolve the rem- 
nant at this wavelength. Due to the combination of increased 
brightness and the ATCA upgrade with the Compact Array 
Broadband Backend (CABB; Wilson et al. 2011), the SNR has 
been resolved for the first time at 3 mm (94 GHz). 

The ability to resolve the remnant at high radio frequencies is 
crucial in investigating the different sources of emission within 
the remnant, and as such the new ATCA observations are a use- 
ful complement to future observations with the Atacama Large 
Millimeter/submillimeter Array (ALMA). 

2. Observations and data reduction 

SNR 1987 A was observed at 3-mm wavelength with the ATCA 
over three 12-hr sessions, on 2011 June 30, July 1 and August 
2 (programme C2495; RI. J. van Loon) - days 8892, 8893 and 
8925 since explosion - with the five antennae equipped with re- 
ceivers operating in the 3-mm band (85-105 GHz). The array 
was in the 7 5 OB array configuration on June 30 and July 1, with 
baselines between 61 and 765 m, and in the HI 68 array con- 
figuration on August 2, with baselines between 61 and 192 m. 
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Fig. 1. DifFraction-limited Stokes-I continuum image of 
SNR 1987 A at 3 mm from observations made on 201 1 June 30, 
July 1 and August 2. The image is restored to a Of'7 circular 
beam (plotted in the lower left comer). The off-source r.m.s. 
noise is ^ 0.5 mJy beam"^. 



Observations were performed in 2-GHz wide bands centred on 
frequencies of 93 and 95 GHz. All three sessions were charac- 
terised by stable atmospheric conditions, with very low precip- 
itable water vapour on June 30 and July 1. 

The bandpass calibrator, PKS B 192 1-293 was used on June 
30 and August 2, while PKS B0637-752 was used on July 1 . The 
phase calibrator PKS B0530-727 was observed for 1.5 min per 
2 min integration time on the source, while the pointing calibra- 
tor PKS B0637-752 was observed for 2 min at approximately 
hourly intervals. Uranus was used as flux density calibrator. 
Observations of SNR 1987A were centred on RA 5^35"^27?975, 
Dec -69°16ai:'08 (J2000) as from Potter et al. (2009). 

The data were processed using the miriaeQ package. The task 
ATFix was first used to correct the system temperatures, instru- 
mental phases and baseline lengths. The data were examined, 
and scans during poor atmospheric phase stability were rejected. 
To maximise the signal-to-noise (S/N) ratio, robust weighting 
was used to invert the calibrated uv data at both frequencies, with 
weighting parameter 0.5 (Briggs 1995). A preliminary clean 
model (Hogbom 1974) was constructed by using a small num- 
ber of iterations (viz. 200). The source had suflftcient S/N for 
phase self-calibration, which was performed with a 2-min solu- 
tion interval. Subsequently, deconvolution was performed using 
the maximum entropy method (MEM; Gull & Daniell 1978). For 
the final image the MEM model of the combined datasets was 
restored to a diff'raction-limited circular beam with full- width at 
half-maximum (FWHM) 0'' 7 - i.e. the synthesized beamsize of 
the 750B array - and regridded at a pixel scale of Of'Ol . Figure 1 
shows the Stokes-I continuum image; the off'-source r.m.s. noise 
is ^ 0.5 mJy beam"^ and we determine an integrated flux density 
of SNR 1987 A of 23.7 ± 2.6 mJy. 

To construct a spectral index image, we also reduced obser- 
vations at 3 cm performed on 2011 January 25 (Ng et al., in 
preparation). The reduction procedures included application of 
uniform weighting, one iteration of phase self-calibration over 
a 5 -min solution interval and MEM deconvolution. As for the 
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Fig. 2. IR-radio SED of SNR 1987 A with data from Matsuura et 
al. (2011), Lakicevic et al. (2011, 2012), Zanardo et al. (2010), 
and this work. Power laws are fit to /I > 30 cm, and two diff'erent 
models for dust emission are plotted (cf. Lakicevic et al. 2012). 
The upper limit of 1 mJy is plotted for any emission in the centre 
of the 201 1 image at 3 mm in addition to that extrapolated from 
the synchrotron power-law. 

3-mm image, the 3-cm image was restored with a 0''7 circular 
beam and regridded at a pixel scale of Of'Ol. 

3. Discussion 

Compared to July 2005 (Lakicevic et al. 2011), the flux den- 
sity has increased 2.1 ± 0.4 times (Fig. 2). The new 3-mm mea- 
surement lies on the position expected from extrapolation of a 
power-law fit to data at /I > 30 cm of the form Sy oc y^, with 
spectral index a = -0.7, just as the July 2005 measurement 
at 3 mm (cf. Lakicevic et al. 2012) obeyed a power law with 
a = -0.75 to the contemporary radio data. It is odd that the 7- 
and 12-mm band flux densities are lower than expected from the 
power law obeyed by the longer- wavelength data. Could this re- 
sult from decay of the synchrotron mechanism? For a decay time 
T = ll/(B^y) yr (with magnetic field B and relativistic y factor) 
of ~ 10 yr, and happening around the characteristic frequency 
y = 4 X 10^5/ Hz of ~ 10 GHz, one would derive 5 ~ 0.1 G 
(which is rather high) and y ~ 150 (which is typical). However, 
this scenario would be hard to reconcile with the ever-increasing 
synchrotron flux, and would imply that the 3-mm flux is not due 
to synchrotron emission. 

At 0('7 resolution the 3-mm emission from SNR 19 87 A is 
clearly resolved (Fig. 1), making it the highest radio frequency at 
which this has now been accomplished. The 3-mm image shows 
similarities with the images at 7 and 12 mm (Manchester et al. 
2005; Potter et al. 2009) and at 3 cm (Ng et al. 2008). The emis- 
sion is broadly distributed along the equatorial ring, but unlike 
the optical images the radio appearance is asymmetric with the 
emission peaking in the eastern lobe (Fig. 3). The flux ratio be- 
tween the eastern and western halves of the image is ^ 1.3. By 
fitting an optically thin shell, this asymmetry was found to be 
^ 1.3 at 7 mm (Potter et al. 2009) and ^ 1.4 at 3 cm (Ng et al. 
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Fig. 3. Radial slices through the 3 -mm image at 8 position an- 
gles. Black dots on vertical axes indicate the position of the 
zero for each slice. The top image is overlain with the 70% and 
90% flux-density levels (white contours). The image centre has 
been shifted to the VLBI position of SN 1987 A determined by 
Reynolds et al. (1995) [RA 05^35^27^968, Dec -69° 16^1:^09 
(J2000)]. 



2008). The maximum asymmetry with respect to the VLBI posi- 
tion in the 3-mm image occurs near position angle 90° (i.e. due 
east-west), where it reaches a ratio of ^ 1.6 (Fig. 3). 

The flux density in the centre is ^ 1.6 mJy, but this includes a 
contribution from surrounding emission due to beam smearing. 
The emission peaks at ±0''6 from the centre, at ^ 3 mJy beam"^ 
(Fig. 3). The central intensity level resulting from two identical 
Gaussians of peak intensity 3 mJy and FWHM 0''7 separated by 
1''2 is 0.8 mJy. This leaves room for additional emission in the 
centre also at a level of ^ 0.8 mJy, which may be compared with 
the 2-(T r.m.s. noise of ^ 1 mJy. Given the crudeness of this es- 
timate we settle on an upper limit of 1 mJy at 2-cr significance. 
This is below the upper limit of ~ 2 mJy placed on contributions 
from free-free emission to the mm emission by Lakicevic et al. 
(201 1). Those authors estimated that the cold dust (~ 18 K) seen 




5^35™28.2' 28.0' 27.8' 

Right Ascension (J2000) 

Fig. 4. The 3-mm image is overlain with the contours of the 3- 
cm image derived from observations performed on 201 1 January 
25 (Ng. et al., in preparation), where the contours correspond 
to 20% flux-density intervals. Both the 3-mm and 3-cm images 
have been restored with a 0!'7 circular beam and centred on 
the VLBI position of SN 1987 A determined by Reynolds et al. 
(1995) (blue square). 



at far-IR and sub-mm wavelengths would contribute only ~ 0. 1 
Jy at 3 mm. As the emissivity law is uncertain, Lakicevic et al. 
(2012) considered alternatives, their most extreme model involv- 
ing warmer grains (~ 33 K) behaving like perfect black bodies at 
long wavelengths - e.g., needles or flufl'y aggregates. That model 
need only require < 0.01 M© of dust, and would yield ~ 1 mJy 
at 3 mm (Fig. 2). Thus our current upper limits do not rule out 
such type of dust (and so little of it) to be present. 

The 3-mm image is compared with the 3-cm image in Figure 
4; these form the basis for the spectral index image (Fig. 5). The 
spectral index a was determined from the ratio of the two images 
where the flux levels exceed 3-(T above zero - still, near the outer 
edge of the displayed area the low flux level combined with the 
inevitable imperfections in the image reconstruction can easily 
lead to large variations in the flux ratio between independent im- 
ages. We observe the following: (1) the brightest radio emission 
is characterised by -0.95 < a < -0.75, and (2) the centre has 

a shallower spectrum, with a 0.65, similar to the tentative 

result obtained at cm wavelengths (Potter et al. 2009). 

Difl'usive shock acceleration theory (Jones & Ellison 1991) 
yields a - (\ - y)l2, where y is the power-law index of the 
momentum distribution of the accelerated electrons, N oc p~y 
and depends only on the shock compression ratio, r, as y = 
(r-F 2)/(l - r). Supersonic shocks are easily established in SNRs, 
and lead to 1 < r < 4, i.e. a < -0.5 (r > 4 is possible for 
relativistic and/or non-adiabatic shocks). We conclude that the 
strongest mm-emission, with -0.95 < a < -0.75, arises from 
moderately-weak shocks with 2.6 < r < 3.0. The central region 
could signify the emergence of a pulsar wind nebula with strong 
shocks (r ~ 3.3), or some non- synchro tronic emission - e.g., 

free-free emission with a 0.1 (indeed, the ejecta are clearly 

visible in Hubble Space Telescope images) or emission from 
grains with high emissivity at mm wavelengths (cf. Lakicevic 
et al. 2012; Wickramasinghe & Wickramasinghe 1993), but the 
spectral index value in the image centre is uncertain. 
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Fig. 5. 3 -mm - 3 -cm spectral index image. The spectral index a 
is defined sls Sy oc y^ and was determined from the ratio of the 
3-cm image from observations performed on 2011 January 25 
(Ng. et al., in preparation) and the 3-mm image from observa- 
tions around July 201 1 . Both images have been restored on a 0^'7 
circular beam and centred on the VLBI position of SN 1987 A 
determined by Reynolds et al. (1995) (blue square). The images 
are overlain with contours representing the 70% and 90% flux 
density levels (white: 3 mm; blue: 3 cm). 
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Fig. 6. The 3-mm image of SNR1987A from 2011 (greyscale 
and blue contours; the beam is plotted in the lower left corner), 
overlain with (green) contours of the average of the deconvolved 
0.3-8 keV Chandra X-ray images from October 1999, January 
2000 and December 2000 (Burrows et al. 2000; Park et al. 2002). 
Contours are between 10% and 90% with 10% intervals. 



There is a striking resemblance between the 3-mm image and 
the earliest, deconvolved X-ray images (Fig. 6), taken with the 
Chandra satellite in 1999 and 2000 (Burrows et al. 2000; Park 
et al. 2002; Racusin et al. 2009), when the blast wave had just 
reached the inner rim of the flash-ionized ring. In particular, the 
peaks in the mm emission in the east and south-west correspond 



to similar peaks in the X-ray emission at that time. The X-ray im- 
ages are broad-band, 0.3-8 keV, and were averaged and aligned 
with the 3-mm image; the 1999 image was taken in 0* order 
of the High Energy Transmission Grating, and biased towards 
higher energies than the 2000 images. The correspondence be- 
tween the high-frequency radio structures and X-ray structures 
is remarkable given that radio synchrotron emission depends on 
the magnetic field strength and is highly anisotropic. 

The resemblance between the mm emission and X-rays has 
diminished since 2000, as the X-ray images have become dom- 
inated by soft X-ray emission from hot gas behind the forward 
shock distributed more uniformly along the ring (Racusin et al. 
2009; Zhekov et al. 2010a). The fact that the synchrotron emis- 
sion seen in 201 1 is linked to plasma through which the forward 
shock plowed more than a decade before suggests that particle 
acceleration in SNR 1987 A is sustained at the interface between 
the ionized plasma of the fast wind from the blue supergiant 
progenitor of SN 1987 A, and the denser circumstellar medium 
from the preceding red supergiant stage. A magnetised fast wind 
might also lie at the origin of the bipolar circumburst morphol- 
ogy - indeed, the radio emission extends in the north-south di- 
rection (Fig. 6; cf. the 7-mm images in Potter et al. 2009). 

4. Conclusions 

We have resolved the remnant of SN 1987 A at 3 mm wave- 
length, at a resolution of 0!'7. The image is dominated by syn- 
chrotron emission from just inside the equatorial ring, generated 
by moderately- weak shocks. We note correspondence between 
the 3-mm emission and structures in the hard X-ray emission 
recorded in 1999 and 2000, i.e. more than a decade earlier. These 
observations corroborate the idea proposed by Manchester et al. 
(2005) that the synchrotron emission arises from plasma behind 
the reverse shock (or reflected shock; cf. Zhekov et al. 2010b) 
which moves relatively slowly through relatively hot plasma re- 
sulting in a moderate Mach number and hence a moderate com- 
pression ratio. We set an upper limit on the emission from any 
dust or pulsar wind nebula in the centre, of ^ 1 mJy (2 cr). 
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